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Abstract 

Distributions in jet production often depend on a soft function, 5, which describes hadronic 
radiation between the jets. Near kinematic thresholds S encodes nonperturbative information, 
while far from thresholds S can be computed with an operator product expansion (OPE). We 
design soft functions for jets that serve this dual purpose, reducing to the perturbative result in 
the OPE region and to a consistent model in the nonperturbative region. We use the MS scheme, 
and in both regions S displays the appropriate renormalization group scale dependence. We point 
out that viable soft function models should have a gap associated with the minimum hadronic 
energy deposit. This gap is connected to the leading O(Aqcd) renormalon ambiguity in jet event 
shapes. By defining the gap in a suitable scheme we demonstrate that the leading renormalon can 
be eliminated. This improves the convergence of perturbative results, and also the stability by 
which non-perturbative parameters encode the underlying soft physics. 
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Soft functions play an important role in the study of cross sections close to kinematic 
thresholds, characterized by jets of collimated hadrons with small invariant mass. These 
cross sections are frequently described by factorization theorems involving hard Wilson co- 
efficients, jet functions describing the jets of hadrons, and a soft function S. The hard co- 
efficients and the jet functions are perturbative, while S encodes universal nonperturbative 
information on soft radiation between the jets. The prototype examples are event-shape 
distributions in e + e~ annihilation for large cm. energies Q |TJ, |2], 0], such as the thrust 
T 0, where T = max^ £V |p, • n\/ J2i Wi\ an< ^ the kinematically allowed range is 

1/2 < T < 1. In the threshold "dijet" region of large thrust, T ~ 1, the events are char- 
acterized by two back-to-back jets, and at leading order in 1/Q the factorization theorem 
has two jet functions and one soft function @, HJ. Other examples include distributions 



for jet broadening M, the heavy jet mass ||, and their generalization to angularities [10 



The dijet region also plays a crucial role in event shapes for massive particles, such as the 
invariant mass distribution of jets from top-quarks ||11|| . For applications at hadron colliders 
soft functions which account for initial state radiation are important | 12]| . Finally, for studies 
of weak S-meson decays to jets, soft functions involving the initial state B play a crucial 
role. Examples are B — > X s 7 and B — > X u ev [13, 13, as well as B —>■ X s £ + £~ [^6 , [17 . 
Here phase space cuts enhance the region where the soft function has a large effect. 

Near threshold one can distinguish two regions. Very close to threshold the distribution 
typically shows an enhanced peaked structure, and nonperturbative information in the soft 
function is important for determining the shape and the maximum of the distribution. The 
size of this "peak region" is set by the hadronic scale Aqcd- Next to the peak region the 
distribution typically falls off and shows a tail-behavior but is not yet highly suppressed. 
The dynamics is still dominated by jets and soft radiation, but in this "tail region" the 
leading soft function can be computed perturbatively since it is probed at scales larger than 
Aqcd- In the tail region operators sensitive to nonperturbative physics are power-suppressed. 
Computations of moments involving integrations over both the peak and tail regions can be 
done with this same power expansion. 

Since S encodes different types of physics in the peak and in the tail region, one possibility 
is to make separate predictions for the corresponding cross-sections. However, phenomeno- 
logically it is often desired to treat both regions coherently. In a pioneering analysis of 
e + e~ — > jets [|18j this was handled by implementing a "hard" IR cutoff on the event shape 
variable "e" . Perturbation theory was used above the cutoff and the perturbative corrections 
were frozen below it, with 

R PT (e, Air) =o(e- ^)iff(e) + - e)i^ L (A IR /Q) , (1) 

where Rp^ L contained perturbative results up to two-loop order with next-to-next-to-leading 
log resummation (NLL). The function Rpt{g, Mr) was then convoluted with a normalized 
soft function model S mo d as dictated by the factorization theorem. With a simple choice 
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for S mo d good agreement with LEP data was found for several event shapes. This cutoff 
procedure does not attempt to treat explicitly the renormalization scale dependence in the 
region where the soft function is non-perturbative, nor does it systematically implement the 
perturbative corrections in this peak region. 

The multi-region issue has also been analyzed in the context of 5-meson decays. In 
Ref. |B5) a perturbative tail was glued to the soft-function model, 

S(u>,fi) = S mod (u)) + 9(lj - A - /i/\/e)5'part(^j^) > ( 2 ) 

where S paTt is the "partonic" soft function obtained from perturbation theory and where the 
argument in the ^-function was chosen such that the tail turns on without discontinuity, using 
the condition Spt(u, \/e(u — A)) = 0. This method provides the correct renormalization 
group behavior for the treatment of the tail region at leading order, and is an improvement 
because it allows the perturbative jet function corrections to be incorporated systematically 
in the peak region. Shortfalls are that in the peak region it still hides the dependence on 
the renormalization scale \x in model parameters, and that the perturbative tail is turned 
on by hand at a particular point, rather than allowing it to appear once it dominates the 
non-perturbative corrections. 

In this paper we develop a procedure for constructing soft function models for jets that i) 
reduce to the perturbative result in the OPE region and a consistent model in the nonper- 
turbative region, ii) exhibit the proper renormalization group scale dependence in the MS 
scheme, hi) have a gap associated with the minimum hadronic energy deposit, and iv) are 
stabilized to perturbative corrections by being free from the leading (9(Aqcd) renormalon. 
We show that the soft function gap parameter is essential for removing the renormalon 
ambiguity of the partonic threshold energy order-by-order in perturbation theory. 

Although our procedure is quite general, in order to make all the steps explicit we will 
carry it out in the context of a specific example. We consider event shapes for top-quark 
jets produced in e + e~ — > ti at cm. energies Q 3> m t . The soft function we construct applies 
equally well for massless event shapes in the dijet region, that is, very little of our discussion 
depends on the presence of the top-quark mass or width. We consider the double differential 
top-antitop invariant mass distribution, d 2 a/dM t dM[, where M t j are either in the peak or 
the tail region. In the peak region near the top mass resonance, s t ,t = iMli ~ m t) ~ m t^t 
where T t is the top-quark width, and we have the factorization theorem [JTT| 

d(x pcak tt .^ . (s t -Q£ + \ „ fst- 



dM? dMl 



a Q H(Q, mt, fi) j d£ + d£~B + ( St ^ , //) B_ ( Sf ^ , /i) S np (£ + , £~ /i) , 

(3) 



which is valid at leading order in m t /Q and T t /m t . Here if is a calculable hard coefficient 
and B± are calculable jet functions, whereas S , np (£ ± ,/i) is a nonperturbative soft function 
which peaks for £ ± ~ Aqcd when /i ~ £ ± . In general the convolution probes momenta 
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^ ~ S t,t/Q in the soft function, and large logs in S are avoided by taking // ~ £^ and 
summing large logs in the jet and hard functions. In the peak region s t j ~ Q-^qcd + m t^t, 
so the nonperturbative distribution described by 5 , np (£ ± ,/i) directly effects the differential 
cross section. On the other hand, in the tail region, s t j 3> QAq CD + m t T t , and the dominant 
momenta in the soft function are £ ~ s t f/Q- In the interesting region this is a perturbative 
scale of £ ± ~ 3—30 GeV or larger, depending on the size of Q. The leading order factorization 
theorem in this tail region is 



dM% dNB 



a H(Q, mt, fi) J dt + dVB + ^ ^ , //) £L ( ** ®* , /i) S part (£+, r , //) , 

(4) 



which is valid to leading order in s tt t/Q 2 , rn t /Q, and Aq CD Q / s t) t- Here the partonic soft 
function S , part (£ ± , /i) can be computed as a perturbative series in ot s . Power corrections at 
C^AqcdQ/ s t,t) are determined from S^p in a manner discussed below, while power correc- 
tions at 0(s t j/Q 2 ) involve new factorization theorems containing subleading soft functions 



(which have been worked out for inclusive 5-decays |22|, |23|). 

The soft function carries information on how soft radiation is associated to the definition 
of the invariant mass variables M t j. To be definite we consider hemisphere mass definitions 
where the soft function for both Eqs. (|3|) and @ is Jll 

S^H) = ^Y. 5 ^ + - k t a )^ r -K b m{YnY\Y^^^^ (5) 

c x s 

Here kf a is the total plus-momentum of soft hadrons in X s that are in hemisphere-a, k~ b is 
the total minus momentum for soft hadrons in the other hemisphere. The soft function for 
thrust is related to the hemisphere soft function by 



S T ( T ) = J d£ + de-5(r — J S(£ + ,T) , (6) 

with t = 1 — T, and we emphasize that S(£ + , £~) is independent of the top-mass. In general 
soft functions are matrix elements of Wilson lines, which in our case are 

Y£(x) = P exp [ig J dsn- A s (ns+x)J , Y n \x) = P exp (ig J dsn-A s (ns+x)j . (7) 

In order to predict the invariant mass distribution in the peak and the tail regions we 
would like to connect Eqs. (0) and In this paper we consider the task of constructing 
an appropriate soft-function that contains both S np and S paiTt and which can be applied 
in the peak and the tail region. In order to be useful the result must remain consistent 
for scales \i ~ S t,t/Q, both in the tail region where s t; t S> <?A and in the peak region 
where s t j ~ + QA. We will consider all large logs to have already been summed by 
renormalization group evolution from Q down to these /i's. So the task is to determine the 
soft function matrix element at these /x's, where it should contain no large logs. 
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To begin, consider modeling the soft function by 



s part (£ + -£ + ,r-r,v) s mod (i + ,r) 



(8) 



where 5'p art (£ ± ,/i) is the partonic soft function computed in perturbation theory, and 
>S'mod(^ ± ) is a nonperturbative model function that is /i-independent and contributes only for 
V s1 ~ Aqcd- In Ref. ]17|] an analog to Eq. was used in the study of b — > s£ + £~ to alleviate 
the issues mentioned about Eq. Taking S pauTt to 0(a s ) this formula provided a simple 
way of incorporating the cutoff OPE moment constraints of Ref. |19| in the model for the 
nonperturbative 5-meson soft function. Here we will argue that, suitably refined, Eq. (H) 
can be used to design soft functions for jets that are consistent with the desired properties 
stated earlier. Defining moments 



n,m] 
mod 



d£+dr (£ + nn m s mod (e + ,n 



(9) 



we will demand that S mo d is normalized, S mod 
moments are finite where we have ~ 



n '°' — 1. We will also demand that higher 
(AqcD) n+m for n + m > 0. 
A virtue of Eq. (|j) is that it produces by construction the proper OPE in Eq. (^) when 
used at a perturbative scale \i = fi op ~ s t j/Q ^> Aqcd where £ ± ~ s t) i/Q- To see this recall 
that t- 



Aqcd, and so we can expand S part for £ ± <C to give 



S(£^ ' i Hop) 



S V a,rt{£ ?/ i op)'S' m od 



— 9 Jf* u ) 9 [1 ' 0] 

(]£+ paxt \ ' ^ op > mod 



d 



5'part(^ j H'op) 



[0,1] 
mod 



Q 2 A 2 



QCD \ 



(10) 



Since = 1 we have the desired result that S , (£ ± ,/i op ) = S , part (£ ± , fi op ) at leading power. 
Computing the renormalized soft function in Eq. (|5D to order a s (Fig. [I] with no n /-bubbles) 
it factors as 1 



*Spart J f 1 ) — 'S'part ?/ i )'S'part 5 A 4 ) 



(11) 



with 



7T 



0(1) ln(£//i) 



(12) 



1 We note that the factorized form of the soft function with respect to the two hemisphere light-cone 
variables £ ± in Eq. ([ll]) allows for the possibility to choose two different /i's at which to stop running the 
two jet functions B± in the factorization theorems (||) and (f|). While we do not expect that relation ( |Tl| ) 
is maintained for non-logarithmic corrections beyond the one-loop level, one can prove that the factorized 
form is maintained to all orders as far the scale-dependence is concerned, as in Eq. (15) |24J]. Thus it 
is possible to treat the situation where St and sj are widely separated and to account for the resulting 
non-global logarithms p5| by choosing both renormalization scales differently. 
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We see explicitly that large logs in S part (£ — £, /a) are minimized for fi ~ i — i. Hence when 
I and i are parametrically different it is the larger of the two that is important for the 
proper setting of the renormalization scale in the soft function. This is compatible with 
the expansion in Eq. (pi]). In the convolution with the jet functions in the tail region in 
Eq. (U), the logs in Spart are minimized for /i = /i op , and S , part (£ ± , fi op ) can be determined by 
a truncated series in a s (fi op ). Thus for Eq. (Q) the result in Eq. (||) works at any order in 
perturbation theory. 

We would also like 5'(£ ± , /i) to give a viable model for the peak region S , np (£ ± , /i) in Eq. ([|) 
when it is applied at a low scale /i = fiy ow ^ Aq C d- Here ^ ~ V s1 in 5 , part (£ ± — /t) for the 
convolution in Eq. (Q). This convolution builds the proper /i-dependence into S'(£ ± , /i), since 
the /i-dependence is determined by perturbation theory exactly as in S , part (^ ± , /a). Thus it 
avoids the issue of having a /^-dependence related to the soft function anomalous dimension 
in the model parameters in S mo d- The convolution with 5p art also generates a perturbative 
tail, implying that S'(£ ± ,/i) is not normalizable. To see this define the cutoff moments 



L r L 



s L[n,m]= / d £+ dr (£+) n (r) m S(£ + ,r,fl). (13) 



Using Eq. (||) with Eq. (|12|) and S^®\ = 1 one finds that for L ^> A QCD the normalization 

smo] = 1 + CFa ^ ] {r2 ~ 21n2 © } + • • • ' (14) 

up to terms of 0(a 2 s ) or C(Aq CD /L). Rather than a deficiency, this behavior of S^ 0,0 ! is a 
necessary feature, as it is consistent with the renormalization equations for 5'(£ ± ,/i). Only 
S^od needs to be normalized. 

For the peak region, perturbative improvements to 5 part in Eq. (§) that cause a large 
change to S, could in principle be compensated by changes to the model parameters in S mo d. 
However, it is quite desirable to make S piir t and S mo d as independent as possible, so that the 
interpretation of the model parameters remains unchanged as we perturbatively improve 
Sj>art- A measure for this independence is the convergence of the perturbative expansion 
for Spart at /ti ow . In general the convolution in Eq. (g) generates double logarithmic terms, 
ln 2 (£//i low ) ~ ln 2 (AQCD/^iow) m 5'(^ ± , /^low) where the scale Aq CD is set by parameters in 
'S'mod- The choice of /t] ow should be small enough to avoid these potentially large logarithms, 
but large enough to ensure the validity of the perturbative expansion in a s (fi\ ow ). Thus 
a satisfactory choice of /ai ow might be difficult to find, and requires careful examination. 
To test this issue we can determine the logarithmic series for S' part (^ ± , /i), by finding the 
partonic soft function in renormalization group improved perturbation theory at LL order, 
NLL order, etc. The renormalization group improved S^art satisfies the exact relation 

s part (£ + ,e-,/2) = [d£ ,+ de'~u s (e + -e ,+ ,i2 lf i )u s (r-£'-^ } i2 )s part (£ ,+ ,£ , ~^o), (is) 
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where U s is the LL, NLL, etc. evolution kernel. As indicated this kernel factors in the 
variables £ + and £~ to any order in perturbation theory [p4"| . Using this RG-improved 
Sp a xt(£ ± , fi) the full 5'(£ ± ,/i) in Eq. (Q) also satisfies the evolution equation (|T5D exactly, 
with a /^-independent S m od(-t )• When the logs are small we can expand the RG-improved 
result to a fixed order in a s (/i), and the resulting and S satisfy the RG to this order. 
We will use this truncated version of the NLL series for 5 , part (^ ± , fi) to test for a choice of fi 
which minimizes large logs in the soft function. This will also provide a test for the stability 
of model parameters to the addition of perturbative corrections. 

Lets construct the NLL partonic soft function using a Fourier transform as in |[26|| . At 



NLL order the partonic soft functions factorize S^ 1 (£+,£') = S^ t L (£ + )S^{£-). The 
Fourier transform of S p£LTt (£) = J d£' U s (£ — £' , fi, fi )S paiTt (£' , is a simple product equation 

S pax t(y, fj) = U S {V, yU, yUo) S part (y, fio) , (16) 

where the position space kernel is 

U(y^^ ) = (zy^r { ^ 0) e K ^K (17) 
The LL results for u and K involve rg Usp and (3q and the NLL results involve r™ sp and p\, 



-i-icusp r /-i-icusp n \ i 

£o_ , ( r i Pi \ a s (/i C 

A) 



(18) 



if /x,/x = — -§— -— - r-l-rlnr + - ^ U >- + ^ln 2 r 

/3<J 1<MA0 V r o Po/ 4tt 8tt/3 



which also agrees with Ref. pTfl . Here r = a a (ji)/a a (fio), Cp = 4/3, (3 = 11 — 2/3n/ and 
/?i = 34C^/3 — 10C*a^//3 — 2CVn/ for n/ light flavors, and the one and two-loop terms of the 
cusp-anomalous dimension are T = ACp and r™ hp = 4CV[(67/9 — 7r 2 /3)C J 4 — 10n//9] p8| . 
To obtain a suitable boundary condition to solve Eq. flT5| ) exactly, we note that the series of 
[8(£)\ii k (£/fi)/£] + plus-functions in S , part (^,/i) become a series of ln fc [z y \ie lE \ in S paI t(y, n). 
Thus in position space we can take a boundary condition where all the logs are absent. 
For example, using the LO boundary condition S paTt (y,fi = —ie~ lE /y) = 1 in Eq. flTBj ) we 
obtain S p ^(y,fi) = exp[K(p,, —ie~ lE /y)]. It is straightforward to verify that this partonic 
soft function satisfies the evolution equation in Eq. (|T6|) . Specifying higher order boundary 
conditions for S p3iT t will then properly specify the subleading non-log terms in the series 
for S'part- For instance, S paxt (y,fi = —ie~' yE /y) = 1 — nCpa s (—ie~ lE /yfj,)/8 fixes the NLO 
boundary condition of Eq. (|i~2|). Thus the general solution to Eq. (|i~5|) is 



= \% elVl Wl/, -ie^ E ly) exp [K(ji, -ie~^ /y)] . (19) 

This result allows us to determine the LL and NLL series. Order by order in perturbation 
theory the Fourier transform (FT) can be carried out analytically since 



d k e^ E ( e [Q(f)e~ e]n ( e m~\ 1 



(20) 
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In addition to the leading logs, this inverse Fourier transform gives contributions to non-log 
terms from the expansion of e ElE /T(l — e), which are subleading to the momentum space 
NLL series. As long as such subleading terms are unambiguously defined order by order and 
obey the RGE, one is free to include them in the NLL result. For our purposes we define 
the LL, NLL, etc. results as the resummed series obtained in position space, since it is in 
this space that the evolution equations are the simplest. With the NLO boundary condition 
and NLL evolution we find 



<Spart(^> A*) — 

+ C F (3 



5(1) + 
C 2 



Ti 



-2C 1 



2 48 3 V J ) 1 I 



+ 
C 1 



a 2 s (jj) 



C 2 \2C 3 



3tt 5 



£ 1 +4C 3 £°-^(£)} 



+ 



7T" 



12 

7tt 2 



S3 24 



17tt 2 C3 



■24C 5 )JC 



6 12 
C 3 



+ C 2 P { 

+ Cf (3 2 J- — + -C 1 
M °l 6 48 

+ (c F /3i + ^ rr sp ) { 



7T 

3 ~ '36 
13tt 4 



7vr 2 C3 



£°4 
3 2880 



Ti 



fC 3 7tt 2 



64 



/ 79tt 6 
V20160 



£ 1 + C3 £0 . 



20C; 



5(£) 



71*5(1) i 
3840 J 



(21) 



where D = 1/ li\6(£) In 1 (£ / fx) / (£ / li)~\ . Note that the coefficients for the terms beyond NLL 
order are incomplete, namely a 2 £°, af£ 2,1 ' , and a 2,3 5(£). We show coefficients for these 
terms because of our convention of specifying the series in position space and using the full 
transform to momentum space. To obtain the complete a 2 C° and a 3 C 2 terms we would 
need to include the non-cusp part of the two-loop anomalous dimension. 

Having determined the desired form of S p3iTt in Eq. (||) and a means to test for large logs, 
we now turn to the nonperturbative information in S mo d and the overlap with perturbation 
theory. To satisfy the moment constraints on S^d one can consider a two parameter model 
with exponential tails j|T8"| 



/exp (£ j £ ) 



~ Af(a,b) (£H 



A 2 



A 2 



exp 



-(£ + ) 2 -(£-) 2 -2b£ + t 
A 2 



(22) 



where Af(a, b) ensures f exp is normalized to one, and b ^ controls the noninclusive cor- 
relation between £ + and £~ . Physically the range — 1 < b < is favored WfQ. In the past 



this and other models used for soft functions in jet physics are taken to be nonzero for 
£^ > 0. This is a natural constraint given that it is satisfied to any order in perturbation 
theory for S , part (£ ± ). With > 0, Eq. (|j) enforces £^ > in S'(£ ± , / u). However, a better 
approximation is to take a soft-function with a gap so that the soft-function model vanishes 
for l± < A, 



s mod (I + , t) = f cn M + - a, r - a) 



(23) 
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a) b) c) d) 




FIG. 1: Graphs for the hemisphere soft function with bubble chains. The solid lines denote Y- 
Wilson lines, and the line with ticks is the final state cut which may also cut a quark bubble. 



Here A encodes the minimum hadronic energy deposit in each hemisphere. 2 Since the 
model parameter A ~ Aqcd it has an 0(1) effect in the tail region where the soft function 
is nonperturbative. Among the model parameters A plays a special role because it enables 
a hadronic interpretation for the variables V s1 > A in S' mo( i(^ ± ). 

Through the convolution in Eq. (||) this gap is transferred to give £ > A in S'(£ ± ,/i). 
This transfer relies on the fact that we have a partonic threshold at zero-momentum, i.e. that 
>5part(^ ± — has support only for v 11 > £ ± . However, this transfer is not entirely straightfor- 
ward because in perturbation theory the partonic threshold has a renormalon which yields an 
O (Aqcd) ambiguity in A. In the Borel transform of the hemisphere soft function considered 
here, this renormalon corresponds to a pole at u — 1/2. Since the soft function is universal 
for massless jets and top quark jets this renormalon is also behind the u = 1/2 Borel pole 



identified by Gardi in an analysis of event-shape distributions in full QCD for massless 
partons. The nature of this soft function renormalon is similar to the well known O(Aqcd) 
renormalon of the heavy quark pole mass definition, but is not equivalent to it; rather it is 
specific to the soft function for jets. For example, the u = 1/2 renormalon pole of the soft 
function that occurs in inclusive B decays is solely related to the heavy quark pole mass, 
and is eliminated by switching to a short-distance threshold mass, see for example For 
the case of the top jet event shape distribution considered here the pole mass renormalon 



2 An even more accurate description of the gap would use £ + £ > Tix min < hut here there is a ^ ± beyond 
which S mo d is exponentially suppressed, so the difference to Eq. (p3|) is very small. 
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is contained in the jet functions JTT|, and is of no concern for the construction of the soft 
function. Only gluon fields appear in the matrix element defining our S in Eq. (|5|). 

Using standard renormalon calculus either based on gluon propagators dressed with mass- 
less fermion bubble chains or on the modified gluon propagator 

,5/3- 



e 



q 2 + i0 V V 2 J (-q 2 ~i0) 1+u ' 
the one-gluon exchange graphs in Fig. [I] give the Borel transform 

1\ 8C F e~ 5 / 6 
tt/3o - |) 



(24) 



B 



S^(t + ,r,»)](u*l) = u »(s(£ + )6'(r) + 5'(£ + )5(£-) ) . (25) 



This parameterizes the leading (9(Aqcd) renormalon ambiguity of the tree- level soft function, 
and has the same form as a shift in the zero point of 5(£ + )S(£~) expanded to first order. 
It is also consistent with the result found by Gardi p9[ for thrust, accounting for Eq. (pi). 



Equation (|25| ) can be generalized to soft function diagrams with an arbitrary number of 
gluons with one gluon modified by Eq. (H). Since one can use the soft limit for the modified 
gluon momentum (compared to the momenta of the unmodified gluons) only diagrams where 
the dressed gluon is external need to be considered. The computation of the contributions 
from the dressed gluon then factorizes from the remaining gluons yielding 



B 



W« + ,r,,)](»4) = ^^.(A + A^.r,,). (26) 



This result parameterizes the leading (9(Aqcd) renormalon ambiguity of the soft function 
at any order. The Borel pole at u — 1/2 leads to instabilities in the perturbative predictions 
as we systematically include perturbative corrections to Spart- As we will see below, such 
instabilities are for example reflected in S becoming negative in certain ranges of £^, or in 
an instability of the £ values where S is maximal. Physically, this ambiguity ties together 
the perturbative physics that we aimed to associate with S'part and the hadronic information 
in S mo d, and it must be resolved by experimental information. 

In order to remove the ambiguity and allow for a stable determination from experimental 
data we would like to use a renormalon free scheme for the gap. Thus we take A = A + 5 
where A is a renormalon- free model parameter for the hadronic threshold, and 5 = 61+62+- ■ ■ 
has a perturbative expansion which cancels the renormalon ambiguity in S'part • Shifting 
variables to £^ = £ — 6 we have 



S(£ + ,£~,fi)= I d£ + I df S paTt (£ + -£ + -6,£--£--6,fi)S mod (£ + +6,r+6) 

dt S p ^(£+-£ + -5,£--£--5,ri f exp (£+-A,£--A) . (27) 



00 j —00 

+00 /"+oo 
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To cancel the renormalon ambiguity we must expand Eq. (^7j) in 5 simultaneously with our 
expansion for S pSLlt = S^ Tt + ^ art + . . ., so that 



W^-M = S p ° art (^,/i) + 



){s 2 S° pait (£ ± ^) + 5 1 S 1 part (£ ± ,i,)} 

) ^^part^J AO 



+ 



^part(^ ) Z 1 ) 



+ 



d 2 



d 2 



d d 

dF + de- 
d 2 

f2- 



d£+ 2 dl' 2 d£+d£ 



+ ... . 



(28) 



Here 5{ ~ 0(a l s ) can be defined with any prescription that removes the O(Aqcd) renormalon 
ambiguity, and simultaneously this prescription will define a scheme for the hadronic param- 
eter A. Note that A is renormalization group invariant, thus A inherits a scale-dependence 
if 5 is not renormalization group invariant. Moreover, we note that quadratic and higher 
powers of Si that appear in Eq. ( |2"8"| ) are required to ensure the consistency of the perturba- 
tive scheme. The terms linear in 5i are the ones relevant for removing the leading O(Aqcb) 
ambiguity, having the same form as Eq. (|26|) . 

In order to motivate a definition for a subtraction scheme associated to 5 consider the 
first moment S^ 1 ' * from Eq. (|i~3|) . For now the upper cutoff L is arbitrary. Starting from 
Eq. (f2~7D we use the OPE as in Eq. ((TO) and expand to linear order in 5 to obtain 



S L[l,0] = g 



£[i,o] 

part 



[4M(A) 

L 



q L[l,0] 
^part 







" d 


d ' 


/ dt 


-f d£~ £ + 


1 — oo 


J — oo 


_d£+ 


f d£~_ 



5 



+ d£~ £ + 







+ 







*Spart {£ , 



Sparer,//) 



■ /*) + 5£fl(A) 



_d£+ dt 

where in the second line we dropped a s corrections to the power correction, and here 



(29) 



^rnod(A) 



-oo 

d£ + 



A 



— oo 
+oo 



dt£ + UJ£ + -A,t-A) 



£ fcxp (£ j £ ) 



(30) 



When Spart in the factorization theorem in Eq. (f|) is replaced by the full soft function S, the 
moment S L ^' ^ appears in the small £^ region, and relates the small momentum contribution 
in the leading order factorization theorem with the first power correction. From Eq. (|26|) it 
is clear that there is a C(Aqcd) renormalon ambiguity in S p ^f which should be canceled 
by the 5-term in Eq. (p9|). A suitable form for S to render the leading order factorization 
theorem and the first power correction renormalon free is 



fd£+fd£- £+S pait {£+,£-,fi) 



L L 

Jd£+Jd£~ F 

-oo — oo 



d 



d 



d£+ dt 



(31) 



S paTt (£+,£-,fi) 
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FIG. 2: Soft function S(£ + ,£^fi) as a function of £ = £ + = £~ with \i = 1 GeV, at tree level (solid 
black line), one-loop (dotted red line), one-loop with renormalon subtraction (light solid red line), 
two- loop NLL (dot-dashed blue line), and two- loop NLL with renormalon subtraction (dashed 
blue line). Results are shown for three models: (o, b) = (2.5, —0.8) (left panel), (3.0, —0.5) (middle 
panel) and (3.5, —0.2) (right panel). All models have A = 0.55 GeV and a gap of A = 100 MeV. 



Note that different choices of L correspond to different schemes for renormalon-free gap 
parameters A. Other ways to define 5 are also feasible. From the expression for >Sp ar t given 
in Eq. (|2l|) we obtain 



Si 



2L 



7T 

<X 2 s(fJ>) 



ln£ + l 



7T* 



Cl 



(3q Cf 
2tt 2 



2 L L 48 



-s)ln| + 4C(3) + ^- 



- 12 



cusp 



iln^ + l 
8 L 8 



(32) 



Note that the one-loop 8\ term is exact, while the two- loop term 82 relies on our NLL 
approximation of Eq. fl2~T|). 

Lets examine the impact of renormalon subtractions on the soft function. In Fig. |2| 
S(£ + ,£~,fi) is plotted as a function of £ = £ + = l~ at tree-level (solid black line) and one- 
loop (dotted and lighter solid red lines). Blue dashed and dot-dashed lines are two- loop NLL 
results to be discussed below. We take fi = 1.0 GeV (ct s (/i) = 0.396) and use the soft model 
function of Eq. fl2"2] ) with A = 0.55 GeV, and three different choices (a, b) = (2.5, —0.8) (left 
panel), (3.0, —0.5) (middle panel), and (3.5, —0.2) (right panel). The dotted red line is the 
one- loop corrected soft function prior to renormalon subtractions, with 5i = and A = A. 
The light solid red line is the corresponding result with a renormalon free gap parameter A, 
and subtraction using 5i from Eq. (p2|) . We use L = A as a representative scheme choice, and 
for simplicity have chosen A = 100 MeV. Other values of A simply correspond to a global 
horizontal shift of all curves by the same amount. While the unsubtracted one- loop soft 
functions have unphysical negative values for small £, we see that the renormalon-subtracted 
curves are alway positive. This effect of the renormalon subtraction is very general, we have 
checked that it is realized for any choice of model parameters, renormalization scale fi, and 
scheme parameter L > A. We illustrate this in Fig. by showing soft functions S(£, £, jj) with 
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FIG. 3: Dependence of the soft function S(£, £, fi) on the renormalization scale [i and the renormalon 
subtraction scheme-parameter L for the model with A = 0.55 GeV and (a, b) = (3.0, —0.5). Lines 
use the same conventions as for Fig. [2|. As indicated the upper and lower panels represent curves 
for fj, = 1.0 and 1.3 GeV, while the left, middle and right panels refer to L/A = 0.5, 1.0 and 1.5. 



A = 0.55 GeV and (a, b) = (3, —0.5), for different choices of [i and L. For the upper (lower) 
panels [i = 1.0 (1.3) GeV, and for the left, middle and right panels we have L/A = 0.5, 1.0 
and 1.5. Note that the soft function has an anomalous dimension, see Eq. ( [To] ) and (|T^), so 
its shape and normalization change when varying //. 

In Fig. the subtracted curves also show a somewhat smaller correction to the £ value 
where their maximum is located than the unsubtracted curves, but this effect is more depen- 
dent on the choice of parameters, such as the L value, see Fig. || At 0(a s ) the perturbative 
series for the peak position has not yet approached its asymptotic behavior, but we expect 
the improvement in convergence for the peak position of the soft function to become more 
pronounced when higher order perturbative results for the soft function are considered. 

To test whether S part suffers from large logs for particular values of fi, the 0(a1) NLL 
predictions for the soft function from Eq. (^l|) are shown as the blue dot-dashed and dashed 
lines in Figs. ^| and |3|. The dot-dashed curves do not have renormalon subtractions, and 
again exhibit negative dips. The dashed curve use our renormalon free A, with subtractions 
given by the terms in the last set of square brackets in Eq. (|28D and 5± and 62 from Eq. ([32]) . 
We see that at this order the renormalon subtractions continue to eliminate the negative 
dip at small £ values. The behavior of the peak location for the two-loop NLL result is 
in general not dramatically improved, but this is simply because the 0(ot 2 s ) soft function 
given in Eq. (£H|) is based on a logarithmic approximation in a region where the logs are 
not large, and hence does not contain the large renormalon terms of the full two-loop soft 
function. Finally, for the lower right panel of Fig. ^, we see an indication for an instability 
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M-m t (GeV) M - m t (GeV) M - m t (GeV) 



FIG. 4: Top invariant mass distribution da/dM 2 dM 2 in the peak region as a function of M — mt 
with M = Mt = Mi accounting only for the perturbative corrections arising from the soft function. 
The left, middle and right panel refer to the respective models and renormalon subtraction scheme 
used in Fig. ^ and the same line specifications are employed. 

due to increasing logarithmic terms for fi = 1.3 GeV and L/A = 1.5. For the model function 
of Eq. ( |2"2"D such regions of instability generally arise for larger values of fj,, and increasing 
positive values of b and L/A. This issue might have to be more carefully examined if 
experimental data suggests that such regions of model parameters are favored. 

The impact of the renormalon subtraction is also significant for the differential cross 
section. Let us first consider the peak region based on the factorization theorem (Q). Since 
we only wish to illustrate the impact of the soft function, we use tree-level jet functions 
B ± (s) = l/(s 2 + T 2 t ) for Q/m t = 5, T t = 1.43 GeV, m t = 172 GeV. We also ignore 
common normalization factors, and evolution factors that sum large logarithms down to the 
low renormalization scale fi\ ow of the soft function, since they affect all predictions in the 
same way. (For the case of top-quark jets, a complete analysis including all these terms is 
carried out in Ref. f24| .) Fig. ^ displays this differential cross section for equal invariant 
masses M = M t = Mt over M — mt for the three parameters sets of Fig. ^|. Again we find 
that using a renormalon free gap parameter improves the convergence of the predictions 
and avoids the problem of negative dips in the cross-section. Interestingly, the curves show 
even better convergence compared to the soft function alone, and show nice convergence 
for the peak location. We find that this is true in general and related to the additional 
smearing that is provided by the width of the jet function. These results illustrate that 
the removal of the (9(Aqcd) renormalon contributions in the soft function is essential to 
obtain a renormalon-free mass measurement from the peak position of the invariant mass 
distribution. We emphasize again that the renormalon issue in the soft function treated here 
is entirely independent of the pole mass renormalon problem, which appears in the massive 
jet function and the top quark pole mass. 

Finally, let us examine the tail region of the differential cross section, using again tree-level 
jet functions and equal invariant masses M = M t = Mt and ignoring common normalization 
factors. To be specific we adopt the model with A = 0.55 GeV and (a, b) = (3.0,-0.5). 
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M - m t (GeV) M-m t (GeV) 

FIG. 5: Top invariant mass distribution da/dM 2 dM 2 in the tail region as a function of M — nit, 
with M = M t = Mf, jU = (M 2 - mf)/Q, and A = 100 MeV. In a) results are shown without 
renormalon subtraction (5 = 0), and in b) with a renormalon free gap parameter A. In a),b) 
we show: three curves using S pa , r t at tree, 1-loop, 0(a 2 ) NLL (long dot-dashed black, medium 
dot-dashed red, short dot-dashed blue), and three curvies using the full S at tree, 1-loop, 0{a 2 s ) 
NLL (solid black, long dashed red, short dashed blue). The latter three curves use the model 
with A = 0.55 GeV and (a,b) = (3.0,-0.5), and reflect the effects from power corrections when 
compared to the former three. In c) and d) we show the difference between the first and second 
set of three curves from a) and b) respectively, at tree (solid black), one- loop (long-dashed red), 
and 0{a 2 s ) NLL (short-dashed blue). 

In Fig. |5|a the tree-level (black lines), one-loop (red lines) and two- loop (blue lines) cross 
sections are shown without renormalon subtractions as a function of M — m t . We use 
/i = (M 2 — m 2 ) /Q to avoid large logs in the soft function when plotting over a wide range of 
scales. The dot-dashed lines use the leading order result in Eq. (f|) with only the partonic soft 
function and no gap, and the solid and dashed lines use the full soft function S from Eq. (g) 
instead and take A = 100 MeV. For a given order in a s the difference between the curves in 
Fig. [5|a reflect the typical size of power corrections, and are plotted in Fig. |5|c. In Fig |5|b the 
same tail distributions as Fig. |^a are displayed, but now with the renormalon subtraction. 
Since the perturbative contributions in .S^t are at the scale fi op it is mandatory to choose 
L of order fi op to avoid large logarithmic terms, as can be also seen from Eq. (|3lD , and we 
adopt the specific scheme choice L = ji op . Comparing the curves in Figs. |5]a,b we see that 
the renormalon subtraction substantially improves the perturbative convergence. Figure [|d 
shows the difference between the solid/dashed and the dot-dashed curves from Fig. ||b. 
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Comparing it to Fig. |5|c we see that the renormalon subtractions lead, as anticipated, to a 
significantly better perturbative behavior for values one would extract from the data for the 
power correction. 3 This illustrates that the renormalon subtracted predictions are essential 
for extracting stable and renormalon-free model parameters from experimental data. A 
scheme such as the one used here, where L = fi, works well for both the tail and peak 
regions, avoiding large logs. If a result for the gap model parameter is determined from data 
in a scheme where L = /i, then Eq. (^) can be used to relate the result to other schemes, 
such as for L = ji/2. 

To conclude, we have provided a prescription for designing soft function models in jet 
production, that can be applied both in the peak region where the soft function is non- 
perturbative and in the tail region where the soft function can be expanded with an OPE. 
The method entails the convolution of the partonic soft function with a normalized model 
function that encodes the nonperturabive information, Eq. (|8]) . It automatically implements 
consistent renormalization scaling behavior in the MS scheme, making the design particu- 
larly useful when dimensional regularization is employed for perturbative calculations. As a 
novel feature we argue that the soft function models need to exhibit a gap which accounts for 
the fact that for real hadrons there is a minimal hadronic energy. This gap is also required to 
devise a systematic scheme to remove the leading O(Aqcd) renormalon that is contained in 
the partonic soft function. In Eqs. (p7| , |28| , pTj ) we have provided a simple definition for such 



a scheme and demonstrated that the removal of the renormalon avoids large uncertainties 
in predictions of the soft function and hence the cross-section in the peak region. In the tail 
region it also reduces the size of fluctuations in the power corrections, since they are other- 
wise affected by the C(Aq CD ) renormalon. It is possible to generalize our method to treat 
also sub leading (9(Aq CD ) renormalons with n > 1, which are expected to have smaller effect 
on the soft function stability. Subtraction of these subleading renormalons might improve 
the numerical stability at higher order in perturbation theory of model parameters in S^od 
not related to the gap. 
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3 Note that our choice of a gap of A = 100 MeV shifts all curves in Fig. ||a,b that use S to larger values of 
AI — m t . This is a significant power correction, it increases these cross-sections by ~30%. However, the 
choice of A does not effect the impact of the renormalon subtraction. 
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